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.We provide the results from a spectral analysis of nuclear decay data displaying annually varying periodic fluctuations. The analyzed 

1—1 data were obtained from three distinct data sets: ^^Si and ^^Cl decays reported by an experiment performed at the Brookhaven 

^ National Laboratory (BNL), ^''Mn decay reported by the Children's Nutrition Research Center (CNRC), but also performed at 

I BNL, and ^^^Ra decay reported by an experiment performed at the Physikalisch-Technische-Bundesanstalt (PTB) in Germany. All 

' Q three data sets exhibit the same primary frequency mode consisting of an annual period. Additional spectral comparisons of the data 

^ to local ambient temperature, atmospheric pressure, relative humidity, Earth-Sun distance, and their reciprocals were performed. 

^ No common phases were found between the factors investigated and those exhibited by the nuclear decay data. This suggests that 

teither a combination of factors was responsible, or that, if it was a single factor, its effects on the decay rate experiments are not 

^—1 a direct synchronous modulation. We conclude that the annual periodicity in these data sets is a real effect, but that further study 

t> involving additional carefully controlled experiments will be needed to establish its origin. 
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1. Introduction 

It is generally believed that nuclear decay rates are constant 
under virtually all conditions uJ] and follow the familiar expo- 
nential decay law. 



N(t) = 



dN_ 
'dt 



-ANoe 



(1) 



where N{t) is the number of atoms remaining at time t, No = 
A^(0) is the initial number of atoms, and A = ln(2)/ri/2 is a con- 
stant determined by the nuclear half- life Ti/2- However, recent 
work has uncovered evidence of an unexpected periodicity in 
several nuclear decays 1 2, 3^ 4^, ^ ^ . An examination of the pe- 
riodic behavior exhibited by the data indicated that the primary 
frequency observed in each data set corresponded to a period of 
close to one year [7, 8J. This observation naturally raises the 
possibility that either the detectors or the decays in question are 
being affected in some fashion by an external influence such as 
seasonal temperature variations or solar radiation. 

The reported periodic variations in decay rates motivated 
other experiments and analyses, which both supported and 



questioned the original results H [lOl [ill \lM- Data obtained 
during the solar flare of 2006 December 13 exhibited a signifi- 
cant dip in the counting rate of ^'*Mn nearly coincident in time 
with a solar flare, thus supporting the suggestion of a connec- 
tion between nuclear decay rates and solar radiation [9]. On 
the other hand, in an analysis of data from several decay mea- 
surements, Norman et al. lllOll reported no evidence for a devi- 
ation from Eq. ([TJ. Neither did Cooper fll], in an analysis of 
data obtained from the radioisotope thermoelectric generators 
(RTGs) employed in the Cassini mission. More recently a de- 
tailed analysis of the response of the detector systems employed 
in Refs. [21] and ||4|] concludes that the annual fluctuations in de- 
cay rates reported in these references are not likely attributable 
to environmental influences, such as seasonal temperature and 
pressure variations thereby suggesting the possibility the de- 
cays, and not the detectors, are affected by an external influ- 
ence II12I1 . 

The objective of the present paper is to present a more de- 
tailed spectral analysis of the decay variations reported earlier. 
Here we study both the decay rate frequency content and, for 
the first time, the phases of the observed signals. In addition to 
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verifying the annual periodicity by three different analyses, we 
show that these phases provide an important tool in discrimi- 
nating among possible explanations of the observed annual pe- 
riodicities. As we discuss below, the new phase information has 
two important implications: phase can have a significant im- 
pact on the Pearson correlation coefficient, r, and chi-squared, 
Xy, values used by Norman et al. 1 10], and properly incorporat- 
ing phases can influence the conclusions to be drawn from their 
data. 

The first experiment analyzed was the determination of the 
half-life of ^^Si, a y6 -decay, performed at Brookhaven Na- 
tional Laboratory (BNL) in Upton, New York, between 1982 
and 1986 L2J- In the BNL experiment, it was necessary for the 
counting rate to be directly measured over several years due to 
the long half-life of ^^^Si, which was determined by the BNL 
group to be Ti/a - 172 ± 4 years. The "'^Si sample had a cop- 
per foil which strongly attenuated the relatively low-energy fi- 
emissions from the parent decay. For this reason, the subse- 
quent decay of the higher-energy daughter species, ^^P (a jT- 
decay with a half-life of 14.3 days), was measured and (assum- 
ing secular equilibrium) was used to infer the decay rate of the 
^^Si. Although insignificant for the power spectrum analysis, 
this decay sequence may play a role in the overall phase of 
the observed oscillations if the two nuclei are affected differ- 
ently. As is typical for experiments of this nature, the nuclear 
decay of interest (in this case ^^^Si) was continuously monitored 
in parallel with a much longer half-life species (^^Cl) whose 
decay rate was assumed constant because of its long half-life, 
7^1/2 = 301,000 years. Both nuclides were measured using a 
gas proportional detector, which alternated between the two nu- 
clides by use of a precision sample changer ||2|,[l3|]. The BNL 
team noted that the ratio of the two counting rates significantly 
reduced apparatus-dependent systematic effects such as drifts 
in the measuring electronics, or environmental changes such as 
temperature, atmospheric pressure and relative humidity. Peri- 
odic annual deviations were observed in the ^^Si, ''^Cl count- 
ing rates, and interestingly in the ratio ^^Si/^^^Cl. The analy- 
sis of possible laboratory temperature and humidity effects led 
them to conclude the deviations were of unknown origin. We 
were able to obtain from the BNL group the original raw data 
and laboratory notes for the entire experiment il4|] . Although 
the experiment ran continuously, measurements were made at 
semi-regular intervals that ranged from 10 hours to 53 days 
with a median interval of 1 day. This non-uniform sample rate 
played an important role in the spectral analysis of the data. 
See Figure[T]for the normalized, undecayed ^^^Si/^*C1 counting 
rate. A detailed discussion of the BNL experiment is given in 
Ref. H. 

The second experiment analyzed the caUbration of a ^^^PuBe 
alloy neutron total body irradiator ISy. Since these data origi- 
nate from an experiment utilizing the decay of ^^^Pu, they are 
directly relevant to the Cooper analysis of ^^'^Pu powered RTGs. 
Although the experiment was conducted (from 1977 to 1987) at 
BNL, it was overseen by Ellis \3\ affiliated with the Children's 
Nutrition Research Center (CNRC) in Houston, Texas. In or- 
der to quantify the neutron flux of the ^^^PuBe sources used for 
the neutron activations, the experiment utilized the activation 



of ^^Mn to ^*Mn, and compared the ^*Mn decays (by observ- 
ing the 846.8 keV y-ray emitted by ^^Fe) to a ^^^Cs standard. 
Here, the unexplained fluctuations were observed in the ^''Mn 
decay data (a yS -decay) and exhibited a seasonal difference of 
approximately 0.5% between the winter and summer months. 
In addition to the consideration of changes in counting efficien- 
cies, and the contribution from the low energy component of 
the neutron spectrum, all other possible causes of the seasonal 
variation focused on effects that would impact neutron emis- 
sion rates. Ultimately these were excluded since none could 
produce the observed variation, nor explain why no oscillations 
were observed in ^Cs. Since no known explanation could be 
provided for the periodic signal, a discussion of the medical im- 
pacts of subjects entering into treatment at different times of the 
year was made. As was the case for BNL, the experiment ran 
continuously, but calibration measurements were made at semi- 
regular intervals that ranged from 5 to 829 days, with a me- 
dian interval of 7 days, and a sizeable gap from 1980 to 1982. 
See Figure[T]for the normalized, undecayed ^^Mn counting rate. 
Details of the CNRC experiment are given in Ref. ||3|] . 



The third experiment dealt with the long-term stability of 
nuclear decay detectors, along with the determination of the 
half-life of europium, and was performed at the Physikalisch- 
Technische-Bundesanstalt (PTB) in Brunswick, Germany, be- 
tween 1983 and 1998 fl. Although die raw '^^Eu, '^'^Eu, and 
'^^Eu data were not available, we were able to obtain data for 
the ^^^Ra used as the long-lived comparison standard {T\i2 = 
1600 yr), and these data exhibited annual periodic fluctuations 
similar to those seen in the BNL and CNRC data. For the ion- 
ization chamber investigations, the experimental apparatus al- 
ternated among the europium samples, ^^^Ra, and no source, 
in turn. This methodology of monitoring the background dur- 
ing each cycle allowed background corrections to be made to 
the data directly. Although ^^^Ra decays via an a-emission, the 
decay chain is much more complicated than in the previously 
discussed experiments, and includes some y6-emissions [15]. In 
addition, the Any high-pressure ionization detector used in the 
experimental setup reported the overall current produced by all 
the emissions in the decay chain. As in the BNL experiment 
above, the oscillation was observed in both the nuclear decays 
of interest (the europium samples) and the long-lived compari- 
son standard (^^^Ra) as well as on both the ionization chamber 
and semiconductor detectors. However, while the oscillation 
persisted in the BNL ratio, forming the ratio '^^Eu/^^^Ra virtu- 
ally eliminated the oscillation. The cause of the oscillation was 
attributed to a discharge effect on detector components caused 
by background radioactivity due to radon and daughter prod- 
ucts which have known seasonal fluctuations. During the 15 
year duration of the experiment, measurements were made at 
intervals that ranged from 1 to 75 days with a median interval 
of 3 days. See Figure [Ilfpr the normalized, undecayed " ^Ra 
counting rate, and Ref. ^ for the details of the PTB experi- 
ment. 



2. Spectral Analysis Methods 



3. Power Spectrum Analysis 



Spectral analysis techniques contain several powerful tools 
for the study of time series, one of which is power spectrum 
analysis 116. 17, 181 . In power spectrum analysis, the power, S , 
at any specified frequency, v (measured in cycles per year), is 
computed from a time series formed from a sequence of mea- 
surements. 

Since all three experiments had non-uniform sampling rates, 
a simple discrete Fourier transform (DFT) could not be utilized. 
To deal with any problems arising from possible time-varying 
nuclear decay rates, three independent spectral analyses were 
carried out (which were led by DJ, PS and RL-i-AL), using dif- 
ferent mathematical techniques. The three different techniques 
employed were the Lomb-Scargle method il9ll20ll . a likelihood 
procedure [21], and Bayesian Spectrum Analysis [22]. In each 
of these methods, the power spectrum analysis results in a Pe- 
riodogram or Power Spectral Density (PSD) which gives the 
power as a function of the frequency for the data. 

One standard technique for addressing whether the data have 
a deterministic component is to estimate the "false-alarm" prob- 
ability i20ll23ll . Using this method, the probability of obtaining 
a power S or larger by chance is given by 



Before applying spectral analysis techniques it is convenient 
to work with the normalized, "undecayed" data, U{t), 



P(false alarm) = 1 - (1 - e"^)*' 



(2) 



where M is the estimated number of independent frequencies 
within a desired bandwidth. For the false alarm probability es- 
timate, the number of independent frequencies is taken to be the 
number of peaks observed in the bandwidth between zero and 
the median Nyquist frequency for each data set. This provides 
a method for identifying the statistically significant peaks in a 
given power spectrum. 

In addition to determining the power as a function of fre- 
quency V, the phase can be determined for a given frequency. 
For V = 1.0 year"', we adopt January 1st as the zero phase ref- 
erence date. It is also important to note that the reported phase 
is an average of the entire data set and that the instantaneous 
phase may be variable. 

When searching for common features in two or more power 
spectra, it is convenient to use the joint power statistic (JPS) 
which is a measure of correlation given by a function of the 
product of power spectra [24]. The JPS has the property that, if 
each power spectrum has an exponential distribution (such that 
the probability of getting the power S or more is e^^), the JPS 
also has an exponential distribution. For the four power spectra 
involved in the experiments, the JPS statistic J is given, to good 
approximation, by I24fl 



3.881X- 
~ 1.269 + X 

where X is the geometrical mean of the four powers. 



X - (S SiS clS MnS Ra) 
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which removes the known decay trend from each data set. In 
the absence of any time-dependence other than that expected 
from Eq. ([T]i, U(t) should be constant. 

All three mathematical techniques were applied to the BNL 
and PTB data, and resulted in the same general conclusion that 
the annual periodicity is a real signal for which the cause is yet 
to be determined. For reference, the PSD's from the Lomb- 
Scargle (LS) spectral analysis of the three experiments are pro- 
vided in Fig. |2] Figure |2] exhibits the peak widths and ampli- 
tudes of the different nuclides and focuses on the annual period. 
The analysis is complicated by aliasing that arises from the in- 
terplay of the very strong periodicity at v = 1.0 year"' with the 
sampling time series. Several independent methods were used 
to identify aliasing peaks in the spectra, and included forming 
the power spectrum of the timing, observing the eff'ects of sub- 
tracting the annual period, and application of reference func- 
tion comparisons. We found that several statistically significant 
peaks at higher frequencies are due to the sample rate and are 
not actually present in the data. The frequency, power, and false 
alarm probability associated with the two most prominent non- 
aliased peaks for each nuclide are given in Table [1] 

4. Phase Analysis 

We can gain significant insight into possible causal factors 
when their expected frequency content and phase information 
are compared to those of the data sets. For example, we would 
expect detector systematic and environmental effects to pro- 
duce the same phase for each sample in experiments in which 
the detector alternated between different samples during each 
run 111211 . From the power spectrum alone, any potential causal 
factor with an annual period would produce a PSD with a 
prominent peak centered at the same frequency as that ob- 
served in the experiments. For possible factors with direct syn- 
chronous modulation, the frequency and phase of the nuclear 
decay data would mirror that of the causal factor thereby pro- 
viding possible rejection criteria. For reference, the phases of 
the four nuclides from the experiments are provided in Table |2] 

Four independent possible factors with known annual period- 
icities and their reciprocals were investigated. These included 
ambient air temperature at the experiment location (T), atmo- 
spheric pressure (P), relative humidity (H), and the Earth-Sun 
distance (R). National Oceanic and Atmospheric Administra- 
tion (NOAA) and United States Air Force (USAF) archived 
observations from John F Kennedy International Airport in 
Queens, New York (near Upton), were used for the BNL and 
CNRC experiments, while observations for Brunswick (Braun- 
schweig) Germany, were used for the PTB comparison. 

Although many possible factors possess an annual period, 
their phases may be significantly different. For example, the 



mid-latitude temperature history of sites in both the northern 
and southern hemispheres would share the frequency content in 
a PSD but would be ~ 180 degrees out of phase. It is impor- 
tant to note that the ambient atmospheric conditions discussed 
here are those of the conditions outside the laboratory and are 
not as stringent as limits that could be obtained if the actual 
laboratory environment had been monitored and reported. For 
the purposes of this paper we make the assumption that these 
represent a conservative estimation, since the laboratory condi- 
tions likely did not experience as dramatic fluctuations as those 
of the outside atmosphere. Additionally, we assume that trends 
within the laboratory generally followed those of the outside 
ambient conditions. For example, we would expect that as the 
absolute humidity of the ambient outside air increased so too 
would the conditions within the lab, albeit to a lesser degree. 
Table [3] provides a summary of the power and phase for each 
factor 

It is natural to search for possible correlations between ex- 
ternal influences and the nuclear decay data using traditional 
parametric tests. We note that the relative phase between two 
periodic functions has a significant impact on goodness of fit 
statistics, including the Pearson correlation coefficient, r, and 
the chi-squared statistic, Xv For example, a periodic signal 
in the nuclear data with a phase different from 1 IR^ could be 
masked in an analysis comparing these data if phase informa- 
tion is not included. This is evident by noting that for two pe- 
riodic signals with a common frequency v, but with phases (pi 
and (^2, '" - cos(0i - 02) when an integral number of cycles is 
compared. For a non-integer number of cycles, the expression 
is more cumbersome but still depends on the relative phase and 
ranges from -1 to 1. For phases on the order of those observed 
in the data the application of r or xl^ as is done by Norman et 
al. [10], would result in a poor fit statistic when compared to an 
otherwise good periodic fit. We show this in detail in Ref. 1251 . 
where we re-examine the ^^Na/**Ti raw data of Ref. jlOll which 
have been generously provided to us. Using the power spectrum 
analysis methods provided above, the ^^Na/^'^Ti data yield a fre- 
quency of 1 .09 + 0.33 year"' with more than a 97% confidence 
the peak is statistically significant and a corresponding phase of 
-25 degrees. When these data are analyzed incorporating the 
30 day phase shift between 1//?^ and the peak counting rates 
noted in Refs. 12,13, Hill, the resulting fit has a helttx xl than 
does the null hypothesis. 

5. Discussion 

We observe from Fig. |2] and Table [1] that all nuclides are 
dominated (to varying degrees) by a strong annual fluctuation. 
This was demonstrated using the Lomb-Scargle (LS) method, 
a likelihood procedure (LIK), and Bayesian Spectrum Analysis 
(BS A), the results of which were in good agreement. For exam- 
ple, the ^^''Ra LS period was 363.8 + 10.1 days, the LIK period 
was 365.3 + 3.3 days, and the BSA period was 363.75 + 0.27 
days. This is also confirmed by the strong power displayed in 
the joint power statistic (JPS) formed from all four nuclei. 

As noted in the Introduction, the phases of the annually vary- 
ing periodic signals may help to discriminate among alterna- 



tive explanations for the observed periodicities. We have de- 
termined the phases of the BNL, CNRC, and PTB data sets 
by three different methods whose results are in close agree- 
ment. Again using the ^^^Ra data as an example, the phases 
of the annual period were -32.0+ 1.6 deg, -30.6 + 3.3 deg, and 
-30.9 +1.1 deg. Representative values using the Nelder-Mead 
(NM) method for this phase are given in Table |2] and in Table 
[3]the phases of possible influences responsible for the periodic 
signals are presented. 

The phases of the ^^Si, ^''Cl, ^^Mn, and ^^''Ra data are differ- 
ent from one another, as well as also different from the possible 
factors T, P, H and R. To illustrate the implications of Tables 
|2]and[3] we note that the phases of ^~Si and ^^Cl are different, 
notwithstanding the fact that both nuclides were observed in the 
same detector in alternating 30 minute time blocks [2]. For ex- 
ample, if temperature variations at the detector were responsible 
for the observed periodic signals, we would expect that the ^^Si 
and "^^Cl phases would have been the same and agree with phase 
of the temperature, contrary to what is observed. By contrast, 
the difference in the observed phases between ^^Si and ^^Cl is 
understandable in a picture in which the nuclides themselves are 
being affected differently. We note that part of the phase differ- 
ence between Si and ''CI may be attributable to the fact that 
what is directly observed in ^^Si are the y6-emissions from the 
decay of daughter nuclei ^^P produced by ^^Si. Since the half- 
life of "'^P is 14.28 days, account must be taken in analyzing 
decay data of a possible phase lag between a change in the ^^Si 
decay rate, and its manifestation via ^ P decay. Specifically, if 
the dominant effect of an external influence on the decay se- 
quence: Si — > P + e + Ve, is on the parent Si nucleus, 
then any modification in the ^^Si decay rate will be reflected 
in the ^^P decay rate (which is what is observed) ~ 20.6 days 
(= T1/2) later This illustrates how a mechanism which focuses 
on the decay process itself (rather than on the behavior of the 
detection system) could account for the different phases that are 
observed in the data. 



6. Conclusions 

Following a spectral analysis of data from three separate nu- 
clear decay experiments (■'^Si and ■"'CI from BNL |2], ^*'Mn 
from CNRC fs], and ^^*'Ra from PTB |4]), we conclude from 
three separate analyses that all exhibit a primary mode with an 
annual period. In addition, we determined the phase of the var- 
ious signals, and found a significant phase shift in each data 
set relative to that from four possible external factors. We have 
shown that ambient outdoor air temperature, atmospheric pres- 
sure, relative humidity, Earth-Sun distance, and their recipro- 
cals could each be excluded as the sole explanation of the ob- 
served annual fluctuations, assuming their effects on the decay 
rate experiments are in direct synchronous modulation. Al- 
though we are unable to use this analysis to prove a given hy- 
pothesis, in light of known environmental detector effects lll2|] . 
and the observation that phase shifts can exist in some exper- 
iments when nuclei are affected differently, we conclude that 
these results are consistent with the hypothesis that nuclear de- 
cay rates may be influenced by some form of solar radiation. 



The inference from our analysis that different nuclei may be 
affected differently by an external source could help to explain 
recent papers by Norman et al. 1 10] and Cooper 1 11] who have 
set limits on possible variations in the decay rates of several nu- 
clides. Since the phases of the various nuclides summarized in 
Table |2]differ sufficiently from that of IjR^, failure to take them 
into account masks the statistically significant annual variation 
present in the data of Norman et al. il0ll25ll . More generally. 



it is reasonable to suppose that the same complex details of nu- 
clear structure (e.g. nuclear wavefunctions, angular momentum 
selection rules, etc.), which are responsible for the fact that half- 
lives vary from fractions of a second to tens of billions of years, 
could also affect the response of different nuclei to some exter- 
nal influence. With respect to the analysis of Cooper Ollll we 
have noted that since the data of Ellis |3] described above also 
derive from a ~^'^Pu source, it is useful to compare these two 
data sets which appear to be in disagreement with each other. 
We observe that in an RTG the decay of Pu produces U 
which has a half-life of 2.46 x 10^ year, and hence an RTG is 
only sensitive to a-decays. In contrast, the ^^^'PuBe source used 
by Ellis was monitored using the y6-decay of ^^Mn. Absent a 
clear understanding of the cause of the oscillations, we sim- 
ply observe that the data of ElUs yfl and Cooper [11] could be 
compatible in a picture in which an external influence predom- 
inantly affects yS-decays. We conclude from the preceding dis- 
cussion that there is no obvious conflict between our results and 
those of Norman et al. 1 10] or Cooper 1 1 1], given our limited 
understanding of the mechanism underlying the annual fluctu- 
ations discussed in the present paper. Additional discussion of 
these papers can be found in Ref. [12]. 

Our analysis can be extended to apply to other possible ex- 
planations for the observed annual fluctuations. As the sen- 
sitivity to nuclear decay rates improves, additional care must 
be taken to monitor room atmospheric conditions, radon levels, 
and power-source fluctuations, as well as the Earth's electric 
and magnetic fields. It is also possible that some mechanisms 
may be excluded by performing experiments where any or all 
are intentionally varied. For example, repeating one of these 
experiments at a mid-southern latitude should flip the phase 
of many atmospheric factors. We thus encourage the develop- 
ment of experiments to reproduce the observed oscillations in 
the specific nuclei analyzed here with an emphasis on carefully 
controlling the external environment. 
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Figure 1: Reproduction of the normalized, undecayed counting rates reported by the Brookhaven National Laboratory fz|. Children's Nutrition Research Center [31, 
and the Physikalisch-Technische-Bundesanstalt [4j. These graphs include the curve fits obtained using the Lomb-Scargle (LS) and Nelder-Mead (NM) methods, as 
discussed in the text. 



Table 1 : Spectral content summary of the BNL, CNRC, and PTB nuclear decay data. Only the two most prominent non-aliased peaks ai'e included for each nuclide, 
and ''^Si/^^Cl denotes the ratio of N(Si)/N(Cl). 



Experiment 


Nuclides 


Frequency [year '] 


Power 


False Alarm Probability 


BNL 


^^Si/'^Cl 


0.98 + 0.12 


49.1 


2.7 X 10""* 


BNL 


3^C1 


0.90 + 0.31 


15.1 


1.9 X 10"^ 


BNL 


32si 


1.10 + 0.23 


15.0 


1.9 X 10"^ 


BNL 


32si 


0.70 + 0.10 


13.0 


1.3 X 10"' 


BNL 


^^Cl 


11.34 + 0.28 


11.0 


9.5 X 10"^ 


BNL 


32si/3^Cl 


0.33 + 0.12 


10.9 


1.1 X 10"' 


CNRC 


^"Mn 


0.98 + 0.04 


17.8 


1.6 X 10"^ 


CNRC 


5^Mn 


0.25 + 0.04 


14.1 


6.4 X lO"'' 


PTB 


^^''Ra 


1.00 + 0.03 


494.9 





PTB 


22^Ra 


0.22 + 0.04 


82.1 





JPS 


All 


1.00 + 0.04 


112.9 





JPS 


All 


0.76 + 0.06 


36.8 


3.7 X 10"'^ 



Table 2: Primary mode phase analysis summary using a 1.00 yeai' frequency for the three experiments. By convention, a phase shift of -41.4 deg coiTesponds to a 
maximum count rate approximately 41 days after January 1. The indicated phases represent the average over the entire data set. 



Experiment 


Nuclide 


Phase [deg] 


BNL 


^^Si 


36.4 + 24.1 


BNL 


''c\ 


106.3 + 7.0 


BNL 


^''^Si/^^'Cl 


-41.4 + 2.7 


CNRC 


^^Mn 


16.2+1.1 


PTB 


226Ra 


-32.0+ 1.6 



Table 3: Phase results of the spectral analysis for the primary modes of enviromnental factors by experiment. As in Table |2] a 1.00 year" frequency was used 
and the indicated phases represent the average for the time interval corresponding to the associated nuclear decay data. The phases of the environmental factor 
reciprocals are provided in the parentheses. 



Causal 




BNL 


CNRC 




PTB 


Factor 


Power 


Phase [deg] 


Power Phase [deg] 


Power 


Phase [deg] 


T (l/T) 


602.7 


153.8 (-26.2) 


766.6 155.0 (-25.0) 


1920.0 


160.2 (-19.8) 


P(l/P) 


64.3 


49.0 (-131.0) 


80.2 48.3 (-131.7) 


13.1 


15.6 (-164.4) 


H(l/H) 


34.3 


129.4 (-50.6) 


38.4 132.9 (-47.1) 


761.8 


16.3 (-163.7) 


R{l/R) 


713.3 


176.7 (-3.3) 


1675.2 178.4 (-1.6) 


2730.7 


178.3 (-1.7) 
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Figure 2: Power spectral densities (PSD) of the nuclear decays from the Brookhaven National Laboratory |2], Children's Nutrition Research Center |3I, the 
Physikalisch-Technische-Bundesanstalt |4], and the joint power statistic (JPS) of all the nuclei. The frequency and power of the closest peak to an annual period are 
provided for each nuclide. These PSD graphs were obtained using the Lomb-Scargle method, as discussed in the text. 



